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We calculate the nonequilibrium dynamic evolution of a one-dimensional system of two-component 
fermionic atoms after a strong local quench by using a time-dependent spin-density- functional theory. 
The interaction quench is also considered to see its influence on the spin-charge separation. It is 
shown that the charge velocity is larger than the spin velocity for the system of on-site repulsive 
interaction (Luttinger liquid), and vise versa for the system of on-site attractive interaction (Luther- 
Emery liquid). We find that both the interaction quench and polarization suppress the spin-charge 
separation. 
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I. INTRODUCTION 

While the nonequilibrium dynamic evolution of quan- 
tum systems has long been extensively studiedjA progress 
is hindered by the tremendous difficulties in solving the 
nonequilibrium quantum many-body Schrodinger equa- 
tion. This situation is going to be changed due to the 
progress in experiments and the development in numeri- 
cal methods. 

On the experimental side, the development in manip- 
ulating ultracold atomic gases makes it feasible to study 
strongly correlated systems with time-varying interac- 
tions and external potentials and in out-of-equilibrium 
situations. The high controllability in ultracold atomic- 
gases' systems provides an ideal testbed to observe 
the long-time evolution of strongly correlated quantum 
many-body systems, and to test theoretical predictions, 
such as the Bloch oscillation^ the absence of thermal- 
ization in nearly integrable one-dimensional (ID) Bose 
gases,' 5 and the expansion of BEC in a random disor- 
der after switching off the trapping potential^ These ef- 
forts allow us to study the nonequilibrium dynamics of 
strongly correlated systems from a new perspective. 

Numerically, many techniques have been developed, 
such as, the time-adaptive density-matrix renormaliza- 
tion group (t-DMRG) £ the time-dependent numerical 
renormalization group,— continuous-time Monte Carlo al- 
gorithm^ and time-evolving block decimation method. 8 
Time-dependent spin-density-functional theory (TDS- 
DFT) has been proved to be a powerful numerical tool 
beyond the linear-response regime in studying the inter- 
play between interaction and the time-dependent exter- 
nal potentialities More tests of the performance of TDS- 
DFT will be done in this paper on the polarized system 
with attractive or repulsive interactions. Compared to 
the algorithms, such as the t-DMRG, this technique gives 
numerically inexpensive results for large lattice systems 
and long-time evolution, but with difficulties in calculat- 
ing some properties, such as, the correlation functions. 

The ID bosonic or fermionic systems accessible by the 
present ultracold experiments j 1 1 ' 12 are exactly solvable 



in some cases^ and can be used to obtain a thorough 
understanding of the many-body ground-state and the 
dynamical properties. The nonequilibrium problems in 
ID system are especially remarkable in which the ID 
systems are strongly interacting, weakly dissipative, and 
lack of thermalizationiH The ID systems, belonging to 
the universality class described by the Luttinger-liquid 
theory, have its particularity in its low-energy excita- 
tions, characterized by charged, spinless excitations and 
neutral, spin-carrying collective excitations. Generically, 
the different dynamics is determined by the velocities 
of the charge and spin collective excitations, which has 
been verified experimentally in semiconductor quantum 
wires by Auslaender et al.^ The possibility of studying 
these phenomena experimentally in ID two-component 
cold Fermi gases jia where "spin" and "charge" refer, re- 
spectively, to the density difference and the total atomic 
mass density of the two internal atomic states, was first 
highlighted by Recati et alr^. The different velocities for 
spin and charge in the propagation of wave packets have 
been demonstrated by Kollath et a h 17 ' 18 in a numeri- 
cal t-DMRG study of the ID Fermi-Hubbard model, by 
Kleine et al^ in a similar study of the two-component 
Bose-Hubbard model, and, analytically, by Kecke et al2& 
for interacting fermions in a ID harmonic trap. Exact di- 
agonalization and quantum Monte Carlo simulations are 
also used in studying the spin and charge susceptibili- 
ties of the Hubbard modeh 21 ' 22 Dynamic structure fac- 
tors of the charge density and spin are analyzed for the 
partially spin-polarized ID Hubbard model with strong 
attractive interactions using a time-dependent density- 
matrix renormalization method. 2 "^ The spin-charge sepa- 
ration is well addressed for this system^ We would like 
to mention here that a genuine observation of spin-charge 
separation requires one to explore the single-particle exci- 
tation, which is studied recently in simulating the excita- 
tions created by adding or removing a single particle] 18 : 24 
The nonequilibrium dynamics in ID systems has at- 
tracted a growing attention in the possible equilibrium 
properties after an external perturbation and the changes 
in physical quantities after the quench] 17 ' 25 i 26 The dy- 
namic phase transition and different relaxation behavior 
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are studied with a sudden interaction quencbi^I. The 
relation between the thermalization and the integrability 
in ID system is well addressed.— The real-time evolution 
for the magnetization in the ID spin chain is also studied 
in great details using the t-DMRG. 29 

In this paper, we study the ID system under an instan- 
taneous switching off a strong local potential or on-site 
interactions, namely, a sudden quantum quench is consid- 
ered. The strong local potential creates Gaussian-shaped 
charge and/or spin accumulations at some position in 
space. After the quantum quench, the time-evolution 
of spin and charge densities is then calculated at later 
times. We tackle this problem using TDSDFT based on 
an adiabatic local spin density approximation (ALSDA). 

The contents of the paper are as follows. In Sec. HH we 
introduce the model: a time-dependent lattice Hamil- 
tonian that we use to study spin-charge separation and 
quench dynamics. Then we briefly summarize the self- 
consistent lattice TDSDFT scheme that we use to deal 
with the time-dependent inhomogeneous system. In 
Sec. IIII1 we report and discuss our main numerical re- 
sults. At last, a concluding section summarizes our re- 
sults. 



II. MODEL AND THE METHOD 

We consider a two-component repulsive/attractive 
Fermi gas with Nf atoms loaded in a ID optical lattice 
with N s lattice sites. At time t < 0, a localized spin- and 
charge-density perturbation is created by switching on 
slowly the local potential, such that the system is in the 
ground state of the system with the additional potential. 
At t = + , the localized potential is removed abruptly 
and/or the on-site interaction is switched off instanta- 
neously. This system is modeled by a time-dependent 
Fermi-Hubbard Hamiltonian as follows: 

#(*) = -r£(4-ci+i«r + H.c.) + cr(t)^TWt% 
+ ^v ia {t)h ia . (l) 



Here 7 is the hopping parameter, c\ a (ci a ) creates (anni- 
hilates) a fermion in the ith site (i 6 [l,iV s ]), a =f, 1 is 
a pseudospin-1/2 degree-of-freedom (hyperfine-state la- 
bel), U(t) is the time-dependent on-site Hubbard inter- 
action of negative or attractive nature, and hi a — c} ia Ci a . 
We also introduce for future purposes the local num- 
ber operator hi — and the local spin operator 

The external time-dependent potential Vi a {t) = 
V^ xt 0(— t), which simulates the spin-selective focused 
laser- induced potential. 0(f) is the Heaviside step func- 
tion which relates the quench dynamics to the modifica- 
tion of the local potential. 9(-f) = for t > 0. V£ xt is 



taken to be of the following Gaussian form: 

[i - (N s + l)/2] 2 



Vg* = W^exp 



2a 2 



(2) 



Here W a is the amplitude of the local potential. We 
discuss the system of conserved particle number in the 
canonical ensemble. The number of atoms for spin up and 
spin down is, iV-j- and N±, respectively. The polarization is 
defined as p= (N^—N^/Nf. The on-site interaction and 
W a are scaled in units of 7 as, u — {7/7 an d % = W^/j, 
respectively. 

A powerful theoretical tool to investigate the dynamics 
of many-body systems in the presence of time-dependent 
inhomogeneous external potentials, such as that in Eq. 
©, is TDSDFT^^. based on the Runge-Gross theo- 
rem^ 2 - and on the time-dependent single-particle Kohn- 
Sham equations. The complication of the problem is hid- 
den in the unknown time-dependent exchange and cor- 
relation (xc) potential. Most applications of TDSDFT 
use the simple adiabatic local spin-density approxima- 
tion for the dynamical xc potential ) 30 ' 33 which has often 
been proved to be successful in studying the real-time 
evolution. 31 In this approximation, one assumes that the 
time-dependent xc potential is just the static xc poten- 
tial evaluated at the instantaneous density, where the xc 
potential is local in time and space. The static xc po- 
tential is then treated within the static local spin-density 
approximation. Very attractive features of the ALSDA 
are its extreme simplicity, the ease of implementation, 
and the fact that it is not restricted to mean-field ap- 
proximation and small deviations from the ground-state 
density, i.e., to the linear response regime. The dynam- 
ics induced by the strong local perturbation discussed 
here cannot be dealt with the theory based on the linear 
response while TDSDFT is a good candidate. 

We here employ a lattice version of spin-density- 
functional theory (SDFT) and TDSDFT J£ In short, for 
times t < 0, the spin-resolved site-occupation profiles 
can be calculated by means of a static SDFT. For times 
t > 0, we calculate the time evolution of spin-resolved 
site-occupation profiles nj CT (i < 0) by means of a TDS- 
DFT scheme in which the time-dependent xc potential is 
determined exactly at the ALSDA level (details see, Ref. 
[Tolj ) . The performance of this method has been tested 
systematically against accurate t-DMRG simulation data 
for the repulsive Hubbard model. 10 It is found that, the 
simple ALSDA for the time-dependent xc potential is 
surprisingly accurate in describing collective density and 
spin dynamics in strongly correlated ID ultracold Fermi 
gases in a wide range of coupling strengths and spin po- 
larizations. The performance of TDSDFT in describing 
the nonequilibrium behavior of strongly correlated lattice 
models has also been recently addressed in Ref. [9(. 

In this work, we use this method to mainly discuss the 
nature of the interactions on the velocities of the density 
and spin evolution. The spin-charge dynamics after a lo- 
cal quench is discussed in Luttinger liquids (for U > 0, 
gapless spin and charge excitations) and in Luther-Emery 



3 



liquids (for U < 0, gapless charge and gapful spin exci- 
tations). We consider at the same time the influence of 
polarization on the spin-charge dynamics. For attrac- 
tive interactions, we limit our discussion on the weak- 
interaction case because for strong attractive interactions 
we found our SDFT code overestimates the amplitude of 
the bulk atomic density waves, which will greatly influ- 
ence the TDSDFT results based on that. 

Experimentally the strong local potential can be ob- 
tained by a blue- or red-detuned laser beam tightly fo- 
cused perpendicular to the ID atomic wires, which gen- 
erates locally repulsive or attractive potentials for the 
atoms in the wires, corresponding to W a > or W a < 0. 
In this paper, we are interested in the repulsive poten- 
tial for the atoms. The charge and spin densities can be 
observed by using in situ sequential absorption imaging, 
electron beams, or noise interference,— which, in princi- 
ple, gives an unambiguous information on the spin-charge 
separation. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section, we report on the results calcu- 
lated by solving the time-dependent Kohn-Sham equa- 
tions. Mathematically the solution of the time-dependent 
Kohn-Sham equations is an initial value problem. A 
given set of initial orbitals calculated from the static 
Kohn-Sham equations is propagated forward in time. No 
self-consistent iterations are required as in the static case. 

For times t < 0, the system is in the presence of a 
strong local potential, which creates a strong local dis- 
turbance in ultracold gases and makes the total den- 
sity and spin-density distributions in the center of the 
system locally different (up to a few lattice sites). We 
are interested in two kinds of quench dynamics. The 
first one is that, at time t = + , the local potential 
is quenched with the time-independent on-site interac- 
tion U(t) = U . The second is that, at time t = + , 
the local potential is switched off and at the same time 
the on-site interaction is quenched instantaneously with 
U(t) — UQ(—t). After the quench, excitations are pro- 
duced. We concern in this paper the subsequent real- 
time evolution of the spin-resolved densities after the 
quench, n la (t) = (*(t)|n i(r |*(t)) with \^>{t)) the state 
of the system at time t. Charge density and spin den- 
sity are defined accordingly as rii(t) = riif(t) + ni±(t) and 
Si (t) = [n it (t) - n^.(t)]/2. 

If not mentioned otherwise, the numerical results pre- 
sented below correspond to a system with Nf — 30 atoms 
on N s = 100 sites, and with open (hard wall) bound- 
ary conditions imposed at the sites i = and i = 101. 
The external potential is chosen to be spin dependent: 
w-f = — 1 and wi = 0, used to form a local density and 
spin density occupations in the center of the system. 
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FIG. 1: (Color online) Charge rii(t) and spin Si(t) occupa- 
tions as functions of lattice site i and time t for N s = 100, 
— N_i — 15, to-f = —1, w_i = 0, a — 2, and repulsive inter- 
action of it = +2. Top panel: ground-state charge and spin 
occupations for times t < (solid line) and at time t = 5 ft/7 
(dashed-dotted line). Bottom panel: same as in the top panels 
but at time t = 10 ft/7 (solid line) and t — 20 ft/7 (dashed- 
dotted line). The charge and spin densities are plotted in the 
top and bottom of the panel, respectively. The arrows in the 
plot indicate the positions where the wave packets propagate. 
In the inset, we show the velocities of the charge v c (open 
circles) and spin v 3 (solid circles) density wave packets as a 
function of the amplitude of the local potential \wf\. Both 
velocities are increasing functions of |wf|. 



A. u > and p = 

In Fig. [TJ we show results for a spin-unpolarized sys- 
tem [Nf = N\ = 15) with repulsive interaction of u = 2. 
At t < 0, a dominant local charge- and spin-density pro- 
files in the center of the system are generated by the 
strong local potential. After the quench of the local po- 
tential, the charge and spin densities evolve and split into 
two counterpropagating density wave packets. The prop- 
agation in time is in fact due to the nonequilibrium ini- 
tial condition. The charge density evolves with a quicker 
velocity than the spin, which is in agreement with the 
general picture of spin-charge separation^ A qualitative 
analysis based on the continuity equation for the mo- 
mentum density can also well explain the phenomena of 
spin-charge separation^ 

We notice a common feature in almost all the fig- 
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FIG. 2: (Color online) Charge rn(t) and spin Si(t) occupations 
as functions of lattice site i and time t with quenches for the lo- 
cal potential and on-site interaction, i.e., Vi„(t) = V^ xt ®( — t) 
and U(t) = UQ(—t). The other parameters are the same 
as that in Fig. [T] The static density (solid line) is shown 
together with two time shots for t — 5 ft/7 (dash line) and 
t — 10 ft/7 (dashed-dotted line). 



ures in this paper, that is, the spin and charge densi- 
ties have an asymmetric forward-leaning shape. This is 
caused by a nonlinear effect, i.e., the different local ve- 
locities in the center and at the edges. Since the local 
velocity is proportional to the density, the higher den- 
sity in the center gains larger velocity than that at the 
edges, which qualitatively explains why the asymmetric 
forward-leaning shape happens during the density prop- 
agation. For perturbations with small amplitude, the 
charge velocity is studied in details by t-DMRG and 
compared to the Bethe-ansatz results with good agree- 
ment. 1 ' For the strong local potential studied here, the 
spin and charge velocities, determined from the propaga- 
tion of the maximum of the charge and spin wave packets 
away from the center, vary with time. We thus calculate 
and compare the velocities determined at the fixed time 
t = 10 ft/7. In the inset of Fig. [TJ we show the spin and 
charge velocities as a function of the amplitude of the 
local potential \w^\. We find both velocities are increas- 
ing functions of \w^\. For the charge background density 
(~ 0.3) in Fig. [TJ the charge and spin velocities by the 
Bethe-ansatz method are v c — 1.15 and v s — 0.75. In the 
limit of wj- — > 0, but w± = 0, our results give v c = 1.3 
and v s = 0.65. The differences are possibly caused by the 
simultaneous local perturbations in the charge and spin 
densities used here, which break the spin-charge scenario 
and couple the spin and charge modes, similar to the 
effects caused by the finite spin polarization (see Sees. 
III-C and III-D). 

In Fig. [21 we study the local potential quench to- 
gether with an on-site interaction quench, i.e., V% a (f) = 
V^e(-t) and U{t) = UB(-t). We find that, the spin- 
and charge-density wave packets split and counterpropa- 
gate as usual but the phenomena of the spin-charge sep- 
aration completely disappears. That is, the spin and 
charge densities evolve with the same velocity. From 




FIG. 3: (Color online) 3D plots for the Charge density rn(t) 
(Top panel) and spin density Si(t) (Bottom panel) as functions 
of lattice site i and time t (in units of ft/7) for a harmonically 
trapped system with iV s = 200, 7V t = N i = 15, V2/7 = 
5 x 10~ 4 , Wf — —1, Wi = 0, a = 2, and repulsive interaction 
of u = 2. 



the Luttinger-liquid theory based on the bosonization 
method^ or from the Bethe-ansatz solution^ one can 
derive that the spin velocity v s and the charge velocity 
v c satisfy v c = v s = Vp in the noninteracting limit, with 
vp — 27sin(7rn/2) the Fermi velocity. The interaction 
between the different species is one of the important in- 
gredients for the spin-charge separation, which explains 
the suppression of the spin-charge separation after the 
interaction quench. Making use of the techniques from 
the cold atomic gases, two different ways of quenching, 
used in Figs. [JJand[21 respectively, can give a clear sig- 
nal that different collective spin and charge dynamics 
happens when starting from the same initial strong lo- 
cal perturbation. We would like to mention that Kollath 
proposed to repeat the dynamics in Fig. [1] in higher di- 
mensions where no separation of spin and charge should 
be seen^ We notice that in Fig. [5] already at short time, 
some density waves coming from the sharp edges begin 
to influence the charge- and spin-density wave packets 
from the center. At larger time, they will mix with the 
original packets. 

In practice, an additional trapping potential is un- 
avoidable in the present experimental set-ups. We thus 
present our simulations for the system in the presence 
of an additional weak superimposed harmonic trapping 
potential, namely, V^ xt in Eq. @ is changed into, 

u-(«s+i)/2] 2 / N* + 1 \ 2 
V^ = W a e . (3 ) 

Here we take V2/7 = 5 x 10 -4 . The three-dimensional 
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FIG. 4: (Color online) Charge rn(t) and spin Si(i) occupa- 
tions as functions of lattice site i and time t for N s — 100, 
Nf — N± = 15, Wf — — 1, w_i — 0, q = 2, and attractive inter- 
action of it = — 1. Top panel: ground-state charge and spin 
occupations for times t < (solid line) and at time t = 5 ft/7 
(dashed-dotted line). Bottom panel: same as in the top panels 
but at time £ = 10 ft/7 (solid line) and t — 20 ft/7 (dashed- 
dotted line). The inset shows the velocities of the charge v c 
(open circles) and spin v s (solid circles) density wave packets 
as a function of the amplitude of the local potential \wf\. 



(3D) plots of the time evolution of the spin- and charge- 
density wave packets are shown in Fig. [3] From the 
figure, we observe that, in the presence of the harmonic 
potential the charge and spin wave packets are highly in- 
homogeneous, but the counter-propagation and the sep- 
aration of the charge- and spin-density wave packets are 
still visible in the background of the inverted parabola. 



B. u < and p = 

In one-dimensional Hubbard model, away from half 
filling, the spin and charge velocities of the low-energy 
collective excitations satisfy ) 37 ' 38 



v s ,c = v F \ 1 T ■ 

V ttvf 

This gives a qualitative explanation that for the positive- 
ly Hubbard model, the charge velocity is larger than the 
spin velocity, while for the negative-?/ Hubbard model, 
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FIG. 5: (Color online) Contour plots for the Charge density 
rii(t) (Top panel) and spin density Si(t) (Bottom panel) as 
functions of lattice site i and time t (in units of ft/7) for a 
harmonically trapped system with N s = 200, Nf = N± = 15, 
V2/7 = 5 x 10 -4 , Wf = — 1, wi — 0, a = 2, and attractive 
interaction of u = — 1. 



the charge velocity is smaller than the spin velocity. In 
Fig. 21 the quench dynamics for the attractive Hubbard 
model, which belongs to the Luther-Emery universality 
class, illustrates that spin-wave packets evolve with a 
faster speed than the charge branches. In the inset of 
Fig. [U we show the spin and charge velocities evalu- 
ated at t = lOft/7 as a function of the amplitude of the 
local potential \wf\. We notice that an abrupt change 
appears in the charge velocity at \wf\ « 0.55. For attrac- 
tive interactions, Luther-Emery paring induces a promi- 
nent density wave characterized by the dip-hump struc- 
ture. While the charge velocity is determined from the 
propagation of the maximum of the charge wave packets 
located at one of the humps of the density wave. The 
increase in the amplitude of the local potential makes 
the maximum of the charge wave packets move from the 
lattice site i = 44 to i = 39, which explains the discon- 
tinuity of the charge velocity for attractive interactions 
at |u>-|>| ~ 0.55. However, this discontinuous change has 
artifacts because the way of extracting i> CjS used here is 
not an optimum one. In Fig. [SJ we present the contour 
plots of the time evolution of the density and spin pack- 
ets for the system in the presence of a harmonic trapping 
potential with V2/7 = 5x 10~ 4 . The different evolution 
velocities for the charge- and spin-density wave packets 
are clearly visible. 
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FIG. 6: (Color online) The ground-state charge and spin oc- 
cupations as functions of lattice site i and time t for the sys- 
tem of repulsive interaction of u — 2 in the polarized case of 
P = 0.47 (JVf = 22, iV; = 8). Besides the ground-state den- 
sity and spin density (solid line), three time shots are shown 
with t = 5 ft/7 (dash line,) t — 10 ft/7 (dashed-dotted line), 
and t = 15 ft/7 (dotted line). 
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FIG. 8: (Color online) The ground-state charge and spin occu- 
pations as functions of lattice site i and time t for the system 
of attractive interaction of u — — 1 in the polarized case of 
p = 0.47 (Nf = 22, N_i = 8). Besides the ground-state density 
and spin density (solid line), three time shots are shown with 
t = 5 ft/7 (dash line,) t = 10 ft/7 (dashed-dotted line), and 
t = 15 ft/7 (dotted line). 
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FIG. 7: (Color online) Same as Fig. 6 but for the polarized 
system of P = 0.87 (JV t = 28, N± = 2). 



FIG. 9: (Color online) Same as Fig. 8, but for the polarized 
system of p = 0.87 {Nf = 28, N± = 2). 



C. u > and p ^ 

The spin-charge separation in a spin-polarized one- 
dimensional system is quite different from the fully po- 
larized one. The spin-charge-coupled dynamics in a po- 
larized system formulated with the first-quantized path- 
integral formalism and bosonization techniques provides 
us a new non-Tomanaga-Luttinger-liquid universality 
class^ For the Luther-Emery liquid of unpolarized at- 
tractive Fermi gases, the spin and charge degrees of free- 
dom are decoupled. In contrast, in the system with fi- 
nite spin imbalance, spin-charge mixing is found based 
on an effective-field theory for the long-wavelength and 
low-energy properties^ In Figs. |5] and [71 the quench 
dynamics for spin- and charge-density waves is shown for 
the system of repulsive interaction (u = 2) with polar- 
ization of p — 0.47 and 0.87, respectively. For p > 0.47, 
there is only small difference between spin and charge 
velocities. In the case of a large polarization, the same 



propagating velocities for spin and charge are obtained. 

D. u < and p =£ 

The quench dynamics for spin and charge density 
waves of the attractive case for u = —1 is shown in Figs. 
[8] and [9l We find with the increasing of the polarization, 
the spin-charge separation is strongly suppressed due to 
the interplay between charge and spin degrees of freedom. 
Theoretically, for the partially polarized system, the spin 
and charge modes are coupled. In this case, there is no 
strict spin-charge separation scenario, namely, the spin- 
charge separation breaks down. Numerically, we observe 
that, at small polarization the spin and charge wave pack- 
ets still evolve at different velocities although they are 
coupled and influence each other. At large polarization, 
the spin-charge separation disappears and evolves at the 
same velocities for both the repulsive and the attractive 
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systems we studied. 

IV. CONCLUSIONS 

In summary, we have calculated the non-equilibrium 
dynamic evolution of a one-dimensional system of two- 
component fermionic atoms after a strong local quench 
with or without interaction quench by using a time- 
dependent density-functional theory with a suitable 
Bethe-ansatz based adiabatic local spin-density approx- 
imation. A test of the performance of TDSDFT is pro- 
vided for the unpolarized systems with attractive or re- 
pulsive interactions in the presence of a harmonic trap- 
ping potential. Under the same local perturbation, the 
charge velocity is larger than the spin velocity for the 
system of repulsive interaction and vice versa for the at- 
tractive case, which is compatible with the low-energy 
collective dynamics from the Bethe-ansatz solution or 



the bosonization techniques. We found the spin-charge 
separation is strongly suppressed when the interaction 
quench is forced together with the local potential quench. 
Spin-charge mixing is found for the system of polariza- 
tion signaling by the disappearance of the spin-charge 
separation. Numerically we observe that the spin-charge 
separation disappears for large polarizations in both the 
repulsive and the attractive Hubbard model we studied. 
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